Full photon statistics of a light beam transmitted through an optomechanical system 
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In this paper, we study the full statistics of photons transmitted through an optical cavity coupled 
to nanomechanical motion. We analyze the entire temporal evolution of the photon correlations, 
the Fano factor, higher moments of the photon distribution, and the effects of strong laser driving. 
In the regime of single-photon strong coupling, this allows us to predict interesting effects such as a 
transition from antibunching to bunching for larger observation time intervals and photon avalanches 
triggered by multi-photon absorption. 
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Introduction - The field of cavity optomechanics, 
where an optical mode couples to nanomechanical vibra- 
tions, has seen tremendous progress during the past few 
years. It is driven by the goals of probing the quan- 
tum motion of nanomechanical devices, implementing 
ultrasensitive measurements, and exploiting the light- 
mechanics coupling for quantum information processing 
(for reviews see e.g. [TJ [2]). Recently, the mechanical 
degree of freedom was cooled down close to its quan- 
tum ground state using optomechanical sideband cool- 
ing [3[[4]. Thus, it now becomes important to analyze 
probes for the quantum dynamics of these systems. Ex- 
perimentally, the most straightforward access is provided 
by the light field. In this context, the photon statistics 
represents a particularly powerful tool, as it is very sensi- 
tive to interactions and reveals the temporal correlations 
produced by the dynamics. A first study in this direc- 
tion recently led to the prediction of photon blockade [5] 
in optomechanical systems. Here, we will take a signif- 
icant step further by analyzing the full statistics of the 
stream of transmitted photons. In contrast to the earlier 
work, this allows us to study many new aspects, including 
strong laser driving, the full temporal structure of pho- 
ton correlations, the Fano factor, and higher moments 
of the distribution of transmitted photons. All of these 
properties become particularly interesting in the regime 
of single-photon strong light-mechanics coupling, which 
(although challenging) is now being approached experi- 
mentally [3l HJ [SHE]- In this regime, the single photon 
coupling rate becomes comparable to the photon decay 
rate, and a few remarkable phenomena have already been 
predicted. A "classical to quantum crossover" may be ob- 
served in nonlinear optomechanical dynamics [9] , photon 
blockade [5 , non-Gaussian [10] and nonclassical [11] me- 
chanical states might be produced, diverse Schrodinger- 
cat states could be generated [12] [13], and QND photon 
and phonon readout and nonlinearities could be enhanced 
in appropriate two mode setups [14] [15] . 

In this paper, we will first consider the photon statis- 
tics under weak laser drive, to make contact to previ- 
ous works. Using quantum jump trajectory simulations, 
we will obtain the photon correlations at all time scales, 
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Figure 1: (a) The standard optomechanical setup, where pho- 
tons are transmitted and detected by a photon detector, (b) 
Nonlinear level scheme induced by the light-mechanics cou- 
pling. The corresponding eigenstates read |n a ,n&), where n a 
denotes the number of photons and nb the number of phonons. 



discuss the Fano factor obtained from photon counting, 
and point out that there are regimes where the bunch- 
ing/antibunching character depends on the time scale 
considered. Then we will turn to strong driving, where 
multi-photon transitions can become possible. We iden- 
tify a special regime of "photon avalanches" , where pho- 
ton absorption triggers a cascade of photon transmission 
through the optomechanical system. There, the light- 
mechanics interplay is particularly important in affecting 
the dynamics. 

The Model - We consider the standard optomechani- 
cal setup consisting of a laser-driven optical cavity, see 
Fig. [TJl. One of the end-mirrors is attached to a vibrat- 
ing resonator, such that the cavity resonance frequency 
^cav(^) depends on the mirror-position x. The coherent 
dynamics of this setup is described by the standard [I] [2] 
Hamiltonian 

H = -HAa^a-hgo (ft + ft) a) a + htttfb+ha L (a f + a) , 

(1) 

written in the frame rotating at the laser frequency ujl- 
Here, a (b) is the photon (phonon) annihilation opera- 
tor, A = uol — ^cav(O) is the laser detuning from res- 
onance, and go = xzpf ■ (d x oj cav (x))\ x=Q is the single- 
photon optomechanical coupling constant. m, and 
%zpf — yh/2m£t are the frequency, mass, and the zero- 
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point fluctuations of the resonator, respectively, whereas 
(%l is the laser driving amplitude. The full dissipative 
dynamics of the open optomechanical system can be de- 
scribed by a Lindblad master equation 
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where p is the density matrix for the photons and 
phonons. K\ [no) is the photon loss rate through the 
input (output) mirror, whereas Tm is the phonon decay 
rate. The standard Lindblad dissipation superoperator 
reads V[d)p = apa) — \a)a, p/2} , where {a, c} = ac + ca 
is the anticommutator. 

Photons and phonons can be decoupled by applying 
the polaron transformation [I2j [13] H' = U^HU to 



Eq. (1), where U = exp gQa)a(W — b)/Q 



This yields 
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mechanical displacement operator £) = 
b)/f2 now enters the transformed driving 
Note that the elimination of the photon-phonon 
coupling has given rise to a photon interaction term. This 
will lead to an anharmonic photon energy level scheme, 
cf. Fig. [I}}, which drastically influences the intracavity 
photon statistics if > kQ [5]. 

Quantum Jump Trajectory Method - In an experiment, 
one is able to analyze the photon statistics by measuring 
the times t n when a photon leaves the cavity through 
the output mirror. To simulate this detection process, 
we will use the quantum jump trajectory method [T6H2Q] , 
which can be understood as an "unraveling" of the master 
equation. It has recently been employed to discuss single- 
phonon detection in optomechanical systems [21]. 

If a photon was detected (a "photon jump" occurred) 
during a time interval [t,t + 5t], then the information 
about the system is updated according to 
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The probability of this event is given by pj = 
K>o(cL^a)(t)St (Note that we assume unit detection effi- 
ciency) . This corresponds to the unraveling of the second 
line of Eq. 

If no photon leaves the output port, then the state 
p evolves according to the first line of Eq. (J2|, but it 
will have to be renormalized after each time step. This 
evolution describes the increase of our knowledge due 
to the absence of a detection event [18] [22]. Overall, 
we obtain stochastic traces of photon detection events 
which directly correspond to what would be observed 
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Figure 2: Typical quantum jump trajectories for a weakly 
driven optomechanical system. These show the expectations 
for the photon number (a,d), phonon number (b,e), and dis- 
placement (c,f), evolving conditioned on the detection of a 
photon leaving the cavity at a certain time. The concomi- 
tant sudden drop in the radiation pressure force leads to sub- 
sequent mechanical oscillations, relaxing during the damp- 
ing time r^ 1 . (a-c) Sideband-resolved regime (k — H/8), 
with go/n = 4. (b-f) Bad cavity limit (k = 5Q, and 
go/ft = 1/10). [Parameters: blue detuning A = Q — gi/Q, 
coupling go = f2/2, laser drive oll — 5 • 10 -3 n, mechanical 
damping Tm — 10 -3 Q, photon decay m — ko — Size 
of Hilbert space: n™ ax = 4 , n^ ax = 9] 



in an experiment. Thereby, we gain access to the full 
statistics, including arbitrary moments and the complete 
time-dependence of photon correlations. 

Full Counting Statistics - The full counting statistics 
p(N,T$) [23 is given by the probability of measuring 
N photons leaving the cavity in a sampling time interval 
Ts, for a constant rate N = ko n of transmitted photons. 
Here, n is the steady state photon number of the (full) 
master equation (J2|. The fluctuations of the detected 
photon number can be characterized via the counting 
Fano factor 
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(4) where (N m ) = £ 



N=0 



l p(N,T s ). If 7- count (T s ) < 
1 there is antibunching (on that observation time- 
scale), whereas for F coun t(Ts) > 1 the photons are 
bunched. Note that the counting Fano factor T is con- 
nected to the two-photon correlation function g^ 2 \r) = 
(at(0)at(r)a(r)a(0))/n via [23] 



^count(Ts) = 1+7V 



[ TS dr [g^(\r\)-l)(l-\r\/T s ). 

(6) 

The long-time limit of the Fano factor then directly 
yields the zero-frequency shot noise power, S(u = 0) = 
]im Tff -.oo(VarJV)/r 5 = iV\F count (T 5 oo). 
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Figure 3: Fano factor J-'count (Ts) for photon detection events, 
as a function of the sampling time interval Ts. (a) Fano fac- 
tor for small sampling times, where antibunching is observed 
(# (2) (0) < 1 and T{T S ) < 1; cf. Eq. (i)). The black dashed 
line is a plot of the short-time solution of Eq. where (a^a) 
and g^ 2 \0) are extracted from a steady-state solution of the 
master equation |2]). (b) Fano factor for larger sampling times 
Ts- We see that .Fcount (T^ — >> oo) > 1 saturates and indicates 
photon bunching. This is a crucial effect of the mechanical 
motion induced by a photon jump (see main text). [Param- 
eters as in Fig. [2p. Total simulation time QT s i m = 2 • 10 8 , 
number of jumps ~ 

In the following, we will use the Fano factor to analyze 
the photon statistics of an optomechanical system in the 
interesting regime of single-photon strong coupling go ~ 
k (and go ~ O). 

Weak Laser Drive - Here, we start by considering a 
weak drive such that the photon number inside the cavity 
remains small, (eft a) <C 1. This is useful to connect to 
previous results [5], in contrast to which we will be able 
to discuss the photon correlations g^ 2 \r) at arbitrary 
time-delays. 

Fig. [2^i-c shows a typical trajectory in the photon- 
blockade regime [5] , where the intracavity photon number 
{eft a) inside the cavity decreases due to a photon jump 
[22] (it would stay constant for a coherent state). In 
general, this single detected photon must have exerted a 
radiation force up to the time it leaves the cavity. Thus, 
upon detection the expectation for the mechanical dis- 
placement is updated to reflect the value it must have 
had due to this force (sudden jump of (x) in Fig. [2^ to- 
wards a greater displacement). The subsequent oscilla- 
tion decays on the mechanical damping time scale. These 
effects could be observed by pulsed displacement mea- 
surements conditioned on photon detection events. The 
photon number itself quickly increases again due to the 
laser driving and then, more slowly on a scale T^ 1 , settles 
to its limiting value. 

In the following, we will focus on the effects of this 
backaction onto the two-photon correlations for time de- 
lays r / 0. As can be seen in Fig. |2^l, where blue detun- 
ing is chosen, the photon number temporarily exceeds 
its limiting value before settling down. This is because 
the mechanical oscillations bring the cavity closer to res- 
onance. Thus, even though at very short times photon 
antibunching is clearly present, longer observation times 
will feature photon bunching for blue detuning A > 0. It 
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Figure 4: (a) Optomechanical level scheme for a choice of 
red detuning that enables two-photon transitions, (b) Map of 
parameters cxl and go /CI, for which avalanche- like behavior is 
expected (dark blue region). Note that in regions (1), (2) and 
(3) the first, second and third inequality of Eq. (|9| is violated. 
The arrow indicates the trace taken in Fig. [5] [Parameters: 
A - -2go/Q and g^/n = 4] 



should be noted that this is a crucial effect of the mechan- 
ical motion and would not be discovered in any analysis 
of the short-time correlations g^ 2 \r = 0). 

Let us now quantify this by considering the Fano factor 
F C ount(Ts) as a function of the sampling time interval Ts, 
as shown in Fig. [3] For the bad cavity case displayed in 
Fig. [2]i-f, the Fano factor would be close to the value 
expected for a simple coherent laser beam, J-'countloo) « 
1. In contrast, for the strongly coupled optomechanical 
system of Fig. [2^i-c, it shows pronounced antibunching 
at short times (J 7 < 1) and bunching (J 7 > 1) for long 
times, up to Ts —> oo. 

In conclusion, we emphasize that F coun t(Ts) has to be 
analyzed for all sampling times Ts to capture the full 
impact of the optomechanical dynamics on the photon 
statistics. Even the sign of the effect (bunching vs. anti- 
bunching) can depend on the observation time interval. 

Photon Avalanches at Strong Driving - At strong driv- 
ing, multi-photon transitions can take place. There, n 
photons enter the cavity simultaneously, since other pos- 
sibilities are precluded by energy conservation due to 
the optomechanically induced photon nonlinearity. We 
will see that if the parameters are suitable, this can give 
rise to "avalanches" of much more than n photons being 
transmitted. 

Such a situation can be evoked by choosing red de- 
tuning, A = — n • g^/Cl. Then, multi-photon absorption 
of n > 1 photons becomes possible and involves n — 1 
virtual transition steps. To favor this virtual (but res- 
onant) |0) —> \n) transition over the non- virtual (but 
off-resonant) absorption of one photon (|0) — » |1)), we 
have to compare the corresponding transition rates (cf. 

Fig. go 
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Figure 5: Photon statistics in the "photon avalanche regime". 

(a) Fano factor J-" C ount(oo) in the long-time limit, as a func- 
tion of laser drive. The highlighted region is the region where 
multi-photon absorption is favored (cf. Fig. JSJd). At the on- 
set of this regime, the Fano factor J-* C ount suddenly rises. 

(b) Rescaled deviation from Poisson statistics, (c) Normal- 
ized third moment, {{N - (N)f) / (N). (d) Distribution of 
photon counts showing strong deviation from a Poissonian 
distribution with same mean. [Parameters: A = — 2#o/Q, 
go = Q/y/2, k = H/(4>/2), T M = 1CT 3 Q. Maximal size of 
Hilbert space: n™ ax = 5 , n^ ax = 36. (d): a L = 0.14ft, 
QT S = 3000]. 

and 

r|0> ^> —\7 ) K^W (8) 

We impose the inequalities 

r| )^|n) > r| )_^|i> , k < |A n | , and k > r| )_>| n ) (9) 

to favor the resonant transition, to ensure that the inter- 
mediate transitions of the resonant transition are virtual, 
and to prevent heating of the resonator. Upon increas- 
ing the laser drive, at some point n-photon absorption 
will start to be favored over other processes. Neglecting 
the mechanical motion, we expect that these n photons 
will decay out of the cavity on a timescale ~ ft -1 , lead- 
ing to photon bunching. In this case, naive expectation 
holds that .Fcount(oo) « n. Thus, the Fano factor has to 
increase when multi-photon absorption starts to become 
important. 

For the simulations, we will focus on the case n = 2, 
such that n = and n = 2 are resonant, cf. Fig. [4^i, but 
the following explanations should hold in general. Fur- 
thermore, we will focus on single-photon strong coupling 
go > n. The remaining free parameters are the laser 
driving strength o^, and the ratio go/£l. A map of these 
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Figure 6: Typical quantum jump trajectory in a regime where 
avalanches are observable due to two-photon absorption. The 
pattern of vertical lines in (a) indicates the actual photon 
jump events, where bunching is clearly observable. The or- 
ange, dashed line in (c) corresponds to the displacement that 
would make the cavity resonant with the incoming laser. [Pa- 
rameters: as in Fig. [5] and ctL = 0.15 Q. Size of Hilbert space: 
rC ax = 5, nr x = 15]. 

two parameters is shown in Fig. [4]d, where the regime 
fulfilling the inequalities ([9| is highlighted. 

In Fig. [5|i, the Fano factor T coun t(oo) is shown as a 
function of the laser drive. We see that /count ~ 1 for 
weak laser driving (Poisson statistics), while .Fcount(oo) 
increases rapidly beyond the threshold of the laser drive 
predicted by the above analysis. There, multi-photon ab- 
sorption is favored and photons are emitted in bunches, 
cf. Fig. [5]i. In contrast to our naive expectation, 
^ r count(oo) is larger than n = 2. This fact is due to 
the mechanical motion induced by the arrival of multiple 
photons inside the cavity: 

In Fig. [6j we display a typical jump trajectory with 
multi-photon absorption events. Upon detection of a 
photon leaving the cavity, the expected intracavity pho- 
ton number actually increases. This behavior is known 
for situations with photon bunching [22j [24 , and it arises 
because the detection of a photon is most likely due to 
a previous multi-photon absorption. The remaining pho- 
tons are likely to be emitted rapidly afterwards. How- 
ever, the mechanical motion induced by the photon jump 
(originally discussed above in our analysis of the weak- 
coupling case) now can bring the cavity closer into res- 
onance, which allows for the transmission of additional 
photons. In this way, avalanches of more than n photons 
become possible (J r C ount( ( ^ ) ) > n). Note that once a sin- 
gle avalanche process terminates, no photon jumps are 
observed any more and the mechanical displacement is 
able to relax towards its equilibrium value. 

For even larger driving strengths, where the third in- 
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equality of Eq. Q is violated, the average photon and 
phonon numbers as well as .Fcount(oo) increase even fur- 
ther. Additionally, no clear avalanches are observable 
anymore. However, the strong photon noise in that 
case can be explained as being due to the transmis- 
sion modulated stochastically by a fluctuating mirror, 
whose motion has been heated. This can be seen most 
clearly from an inspection of the rescaled Fano factor, 
(^count — 1)M (Fig. [5J3). Note that the statistical data 
obtained from the quantum jump trajectory simulations 
(involving ~ 10 5 photons) also provide access to higher 
moments and the full counting statistics at arbitrary ob- 
servation time intervals. The third moment (Fig. [5^), 
for example, displays features such as a maximum in the 
"avalanche regime" studied here. 

Conclusion - We have analyzed the full statistics of 
photons transmitted through an optomechanical system, 
based on quantum jump trajectory simulations. One 
finds that the photon correlations crucially depend on the 
observation time and may even change from antibunch- 
ing to bunching in the long-time limit. For larger laser 
driving, photon avalanches may be observed, where the 
Fano factor strongly rises, beyond the values predicted in 
a naive analysis. These phenomena rely heavily on the 
light-mechanics coupling and the long relaxation times 
introduced by the mechanical system. 
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